Nu Function

private function Nu(lat, a, e) result(radius)

Compute the radius of curvature of the ellipsoid perpendicular to the meridian at given latitude

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: lat

latitude [radians]
real(kind=float), intent(in) :: a

semimajor axis [m]
real(kind=float), intent(in) :: e

eccentricity

Return Value real(kind=float)


Source Code

FUNCTION Nu &
!
(lat, a, e) &
!
RESULT (radius)

IMPLICIT NONE

!Arguments with intent(in):
REAL (KIND = float), INTENT (IN) :: lat !!latitude [radians]
REAL (KIND = float), INTENT (IN) :: a !!semimajor axis [m]
REAL (KIND = float), INTENT (IN) :: e !!eccentricity

!Local variables:
REAL (KIND = float) :: radius

!------------end of declaration------------------------------------------------

radius = a / (1. - e**2. * SIN(lat)**2) ** 0.5

END FUNCTION Nu